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Abstract. - The ground state susceptibility of a system consisting of three fiux-qubits was 
measured in the complete three dimensional flux space around the common degeneracy point 
of the qubits. The system's Hamiltonian could be completely reconstructed from measure- 
ments made far away from the common degeneracy point. The subsequent measurements 
made around this point show complete agreement with the theoretical predictions which follow 
from this Hamiltonian. The ground state anti-crossings of the system could be read-out directly 
from these measurements. This allows one to determine the ground-state flux diagram, which 
provides the solution for the non-polynomial optimization problem MAXCUT encoded in the 
Hamiltonian of the three-flux-qubit system. Our results show that adiabatic quantum compu- 
tation can be demonstrated with this system provided that the minimal energy gap and/or the 
speed of the read-out is increased. 



In the field of solid state qubits superconducting Josephson junctions qubits are one of 
the most promising candidates for quantum computation [1]. Currently they are attracting 
considerable attention, mostly, because they are potentially scalable, can be accessed rela- 
tively easily and controlled individually. Quantum coherent oscillations and conditional gate 
operation have been demonstrated in a two-qubit system [2,3], entanglement between the 
qubit and detector was achieved recently [4], and a lot of attempts have been made in the 
direction of improving the qubit 's coherent dynamics [5]. These experiments were aimed to 
construct a universal set of gates as a basis for quantum computation. 
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An alternative approach is adiabatic quantum computation (AQC) [6]. It is based on the 
encoding of a non-polynomially hard problem in a complex multi-qubit Hamiltonian. This 
encoding should be done in such a way that the ground state gives the solution of the problem. 
In order to reach this ground state, adiabatic quantum evolution of the systems Hamiltonian 
is used. As a preparation step, the system is moved to a flux configuration where it can 
easily relax to the ground state. In the subsequent calculation step, the system adiabatically 
evolves to the configuration where its ground state encodes a non-polynomial hard problem. 
In the final step, this ground state should be readout. A scalable architecture for AQC using 
superconducting qubits was proposed by Kaminsky et al. [7] and the equivalence of AQC to 
standard quantum computation has been shown by Aharonov et al. [8]. 

For flux qubits, such a controllable Hamiltonian can be built up by making use of the 
flux dependence of their energies. Therefore, the (energy) eigenvalues of this Hamiltonian 
become a N-dimensional function of the flux through each qubit: Ei(f q i, f q 2, / g jv) with 
fqi = &qi/&o — 0.5 being the flux normalized by the flux quantum <f> = h/2e relative to 
the degeneracy point of the qubit. The cigenstates of this Hamiltonian can be described 
as a superposition of the 2 N states in the natural qubit (||) or |f)) basis: W = ax|jij2 
• •• In) + a2|ti|2 ••• In) + ••• + a2 N lTit2 ••• In)- If we set the fluxes far away from the 
degeneracy point of each qubit, i.e. \f qi \ » 0, then the system is neither entangled nor in a 
superposition of states but just in the simple classical state: ctj = 1, a\ = ... = ctj-i = ctj + \ = 
... = a 2 N = 0. If the fluxes through all coupled qubits are close to their degeneracies, then 
we approach the common degeneracy point of the whole system, where its state is maximally 
entangled. Close to this degeneracy point the N-dimensional energy surface has a rich structure 
with significant curvature From this curvature, which is directly proportional to the 

system's susceptibility, the qubits' state can be determined [9]. In this Letter we show the 
complete characterization of the susceptibility of a three flux qubit system around its common 
degeneracy point. 

The susceptibility can be determined by using a radio- frequency tank circuit. In this 
technique the tank circuit plays the role of a parametric transducer, since any variation of 
the real part susceptibility \ 0I ; f° r instance, the superconducting interferometers [10,11] is 
transformed into a change of the tank's resonance frequency. The variation of the resonance 
frequency can be detected by measuring the phase shift between the alternating current J r f 
applied through the coil and alternating voltage over the tank circuit. These quantities are 
related through 

tan0 = -^ X , (1) 
i^x 

where Lx is the tank coil inductance and Qx its quality. 

For an 7V-qubit system in the ground state x — 2 Y],, E r °^E ' where r§ v = 

(j2iLi ^i(0|ci 4 V)) x •\7'H CT ^|0))> = M qiT I p i, M qiT is the mutual inductance 

between qubit i and the tank, I pi is the persistent current of qubit i [9] . The susceptibility x 
is significant near the anti-crossings of the energy levels. These anti-crossings are located close 
to the degeneracies of the classical states (Aj = 0), e.g. -E(TTI) = -E(iti) for three qubits. 
Therefore this method allows one to find the location of classical crossings in the full flux 

( 1 )Usually in the common degeneracy point the spacing between ground Eo and first excited state E\ (energy 
gap) is small in comparison with other flux subspaces. The value of this gap is proportional to the qubits' 
tunnelling amplitudes Aj and inversely proportional to the interaction strengths Ji between them. As the 
excitation probability is inversely proportional to the energy gap the system can get excited at this point 
and AQC will give a wrong answer. Therefore, one should increase the energy gap by optimizing the sample 
parameters. Also the evolution speed should be increased in order to avoid thermal excitations while still 
being slow enough to prevent Landau-Zener transitions. 
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Fig. 1 - (Colour on-line) (a) Electron micrograph of the sample, (b) layout of three coupled qubits. 
Each 14 x 7 fim 2 loop contains six Josephson junctions. Two of them (120 x 2000 nm 2 ) are shared 
and couple the qubit to its neighbors, three junctions (120 x 600 nm 2 , 120 x 390 nm 2 , 120 x 600 nm 2 ) 
form the qubit and one auxiliary junction (120 x 2000 nm 2 ) is used to avoid parasitic junctions. Green 
(light gray) and red (dark gray) layers are bottom and top Al layers respectively. The arrows indicate 
the directions of persistent currents in each qubit loop when they flow trough the coupling junctions 
as indicated by the 0, symbols. Since it is energetically preferable to have the same direction of 
persistent currents through the coupling junction the qubits become coupled in this case. So one can 
check that all qubits pairs are anti-ferromagnetically coupled in this design. 




space of a multi qubit system. For finite temperatures, in addition, the higher states become 
populated thereby allowing the readout of their anti-crossings [12]. 
The effective Hamiltonian for three coupled qubits is 

H=-flh4 ) +^ ) }+ E -vi'M , ( 2 ) 

i=l l<i<J<3 

with ej the energy bias on qubit i, A,; its tunnelling amplitude, Jij the coupling energy between 
qubits i and j, and a z , a x are the Pauli matrices in the natural basis of J.) ((j|er 2 |j,) = —1) and 
IT) ((Tl^zlt) = !)■ Depending on the sign of Jij such a Hamiltonian can describe systems with 
anti-ferromagnetic [13, 14] and ferro-magnetic [15] interactions as well as a system with both 
types of interaction [16]. The sample described here has only anti- ferromagnetic interactions. 
As a result the system becomes frustrated around its common degeneracy point. In fact, 
in the classical limit (Aj — > 0) the Hamiltonian encodes the non-polynomial MAXCFJT 
problem [9,17]. The solution of this optimization problem is the classical state with the lowest 
energy. By using the adiabatic evolution of the quantum system with A,; ^ one can find 
all anti-crossings of the ground state and reconstruct the classical cross-overs from them. As 
was explained in Ref. [9], this information allows one to determine the classical state which 
corresponds to the solution of MAXCFJT problem ( 2 ). 

Figurc^shows the sample measured which consists of three Al persistent current qubits [18] 
placed inside a Nb pancake coil. Two junctions in each qubit are nominally 600 x 120 nm 2 , 
while a third one is ^35% smaller. Each qubit is coupled to the other two both magnetically 
and via shared 120 x 2000 nm 2 junctions [12,19]. In order to prevent the formation of junctions 
in the leads, which would lead to parasitic coupling, we add one auxiliary large junction in 
each qubit circuit. The addition of three large junctions increases the effective inductance of 
the loop and modifies the Josephson potential of the qubits. In order to compensate for this 
effect the size of the smallest junction has been reduced. In comparison with previous designs 

( 2 )as an example, for the present system with f q i = 0.008, f q 2 = 0.02 and f q % = the solution of MAXCUT 
problem is described by the state | tT!)> see fig-Et- 
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Fig. 2 - (Colour on-line) (a) — tan 8(IbT, Ibi) at Jb2 = 400 fiA. The qubits' coupling strength 
Ji 2 // p2 $o = M, 2 t5/ 2 /$o » 0.034 and J 13 // p3 *o = M 53T 5/ 3 /#o « 0.033. (b) -tan(9(I bT ,/ b2 ) at 
Tbi = —400 /iA. The coupling strength between qubits 2 and 3 is J 2 3/-Z P 2 ( I > o = Mg^TSI^/^o ~ 0.034. 



with only one coupling junction [9,12] this allows to couple all pairs anti-ferromagnctically. (sec 
caption of fig. \l\for details). The flux through the qubits can be individually controlled by 
direct currents through the coil ibT, and two additional lines Jbi and I\, 2 . The Nb coil has 
an inductance Lt — 134 nH, and together with an external capacitance Ct = 470 pF forms 
a parallel tank circuit with wt/27t = 20.038 MHz and quality Qt = ojt-RtCt = 700 (here 
i?T is the effective resistance). The qubits were fabricated by e-beam lithography and two- 
angle shadow evaporation, whereas the Nb coil was made by e-beam lithography and CF4 
reactive-ion etching. 

The qubit-coil mutual inductances were extracted from the <J>o periodicity of the ac- 
susceptibility of the individual qubit as M q n « 45.77 pH and M q2 T ~ 46.62 pH, M q3 T ~ 
45.79 pH. Figure [3 plots the tangent of the phase shift — tan#(/bT, Ibj) (j = 1 5 2), as mea- 
sured at the mixing chamber temperature of 10 mK. The Hamiltonian itself was completely 
reconstructed from a number of scans away from the common degeneracy point. Each trace 
corresponds to a single qubit anti-crossing, while the repulsion between them confirms the 
anti-ferromagnetic coupling between all pairs of qubits [12]. The individual qubit parame- 
ters (persistent currents I p i and tunnelling amplitudes Aj) were found from the shape of the 
peaks, — tan0(/bT), when all other qubits were biased away from their degeneracy [20] (as 
an example we can choose the curve — tan0(/bT) for 7bi = 600 and 7b2 = 400 /J.A). The 
coupling amplitude Jy between all pairs of qubits can be obtained by measuring the peak- 
to-peak distance 6 fa = M qi TSIi/$>o [12], see fig. El Thus, the dimensionless coupling energy 
is equal to 6fi, i.e. Jij/I p i$o = Sfi. The reconstruction of the parameters resulted in Ai w 
A 2 « A 3 = 70 mK, J pl « Z p2 = 115 nA, I p3 = 125 nA, J 12 /I p2 <P Q = M q2T 6I 2 /$ Q = 0.034, 
Ji3/I P 3®o = M q3T 6I 3 /<£> = 0.033 and J 23 /I p2 <S> Q = M q2T 5I^/^ = 0.034, at the effective 
temperature T e g = 70 mK of our setup known from previous experiments [12, 16]. The dis- 
crepancy between effective and mixing-chamber temperature is due to the noise from external 
leads and amplifier. The substitution of the persistent currents into the later three expressions 
gives J12 ~ J13 « J 23 — (0.61 ± 0.02) K. The slope of the traces in fig. Ogives the mutual 
inductances between qubits and rfc-bias lines. We found M q \^x — 0.595, Af g2 .bi = 0.387, 
M q3M = 0.329, M ql . h2 = 0.296 and M q2 , h2 = 0.347, M q3lh2 = o'.543 (all in pH). 

The experimental reconstruction of the mutual inductance matrix Af^j bi (bi = bl, b2, bT) 
facilitates the possibility of realising automated flux control for each qubit. Indeed, if we want 
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Fig. 3 - (Colour on-line) Plot of — tan(9(/ ? 2,/ g 3) for different fluxes f q i through qubit 1. The left 
column (a, c, e, g, i) represents the tank circuit phase shift, measured for a mixing chamber temper- 
ature of 10 mK, while the right column (b, d, f, h, j) is the theoretical prediction for a qubit system 
at T c ff = 70 mK. The black dashed lines denote the cross-overs between the different classical states. 
The imperfect compensation visible in the upper left corner is due to the saturation of one of the 
current sources, therefore the points there have a slightly larger |/ q i|. 

to set the fluxes (/gi) /g2> /as) m the qubits, we have to feed the bias lines with the direct 
currents: J°j = $o Yli Mji.bi/gi- Figure |3| shows — tan#(/ g2 , f q 3) for different fluxes through 
qubit 1. The left column (a, c, e, g, i) represents the tank circuit phase shift, measured 
at a base temperature of 10 mK, while the right column (b, d, f, h, j) is the theoretical 
prediction [13] for the sample parameters determined from the fig. |21 The experimental and 
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Fig. 4 - (Colour on-line) — tan8(f q 2, fqi) at f q i = 0.008 for hypotethical qubit parameters: J12 = 
J13 = J23 = 0.3 K, /pi = 130 nA, 7 P 2 = I p s = 180 nA, T e g — 70 mK. All other parameters are the 
same as in experiment. Dashed lines denote the classical cross-overs between the different two-qubit 
states. 



theoretical data are in good agreement. 

Figure |3| can be understood if one looks at the ground state in the classical limit (A x — » 
0, A 2 — * 0,A 3 — > 0). In that case we have sharp crossovers between the basis states. In 
the figures they are indicated by black dashed lines between the states IT1T2 J.3) , ITT}, ITI1) 
and Hit) in fig. and the states ||TI), |IH), |TU) and lilt) hi fig. El- Here the spin-up 
and spin-down notation describes different directions of the persistent current in the qubit 
loops. There is also a crossover from |TTI) an d ITIT) to ||TI) an d lilt) at the plane of fig-EJ 
{fqi = 0). The ferromagnetic states ||TT) an d 1141) arc n °t reached for this flux subspace, 
because of their relatively high energies. 

For finite A, macroscopic quantum tunneling removes the classical degeneracy and the 
qubit 's wave function becomes a superposition of the spin states in the vicinity of the anti- 
crossing. As a result the ground and first excited state of the multi-qubit system exhibit 
significant curvature. In accordance with formula the phase shift is maximal at the anti- 
crossings and from the measurements one would be able to reconstruct the classical anti- 
crossing curves. Nevertheless, this is not the case for all parameter space because we measure 
the susceptibility with respect to the total flux. For the transitions from ||TT) to |Ttl) or |Tlt) 
the total "magnetization" does not change, thus there is no susceptibility change and no phase 
shift. Moreover, for our set of parameters, in the vicinity of these transitions the difference 
between the ground and first excited state is smaller than the system's effective temperature. 
Therefore the ac-susceptibility is suppressed because of the partial occupation the first excited 
state [12]. Nevertheless, the transitions would be visible if the persistent currents were different 
as shown in fig. In this case the magnetization of the states mentioned above would be 
slightly different resulting in a change of the susceptibility at their crossover. 

The very good agreement between the measured susceptibility and its prediction underlines 
the usefulness of the parametric transducer for characterising qubits in the ground state. By 
using this equilibrium measurement technique we can find some of the anti-crossings for the 
flux space | f q \ |> 0.003 and all for | f q i |< 0.003. The range where complete characterization 
is possible can be even increased by either decreasing the effective temperature of the setup, 
or by fabricating another sample with a higher A/ J and qubits with different critical currents. 
As an illustration we calculated the system response (see fig. 0} for hypothetical parameters 
J = 0.3 K, ipi = 130 nA, I P 2 = I P 3 = 180 nA and with all the others the same as for the 
current sample. Here all classical crossovers are clearly visible even for T c g = 70 mK. By using 
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non-equilibrium measurements with adiabatic evolution it should be possible to determine the 
state completely. In this kind of measurement one would prepare the system in the ground 
state at a point where this state is easily reached and, subsequently let it evolve to the problem 
Hamiltonian adiabatically slowly but fast enough to avoid thermal excitation ( ). 

In conclusion, the three flux qubit susceptibility was completely reconstructed. The ex- 
perimental data are found to be in complete agreement with quantum mechanical predictions 
in full parameter space. We also demonstrate a controllable multi-qubit ground band anti- 
crossing read-out, which allows ground-state computation with superconducting flux qubits. 
The next steps in this direction would be the improvement of the read-out speed and demon- 
stration of adiabatic quantum computation and its efficiency. 
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